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Abstract 

An improved composite-boson theory of quantum Hall ferromagnets is formulated 
both for the monolayer and bilayer systems. In this scheme the field operator describes 
solely the physical degrees of freedom representing the deviation from the ground state. 
Skyrmions are charged excitations confined to the lowest Landau level. By evaluating 
the excitation energy of one skyrmion in the interlayer-coherent phase it is shown that 
the bilayer QH state becomes stabler as the interlayer density difference becomes larger. 

PACS: 73.40.Hm, 73.20.Dx, 73.40.-c, 75.10.-b 
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1 Introduction 

The quantum Hall (QH) effect is a remarkable macroscopic quantum phenomenon in the 
two-dimensional electron system [E. The underlying physics is understood by the composite- 
boson (CB) picture ^ f|, Q or by the composite- fermion picture || 0. There are several 
new approaches || || to the problem by extending these pictures. 

When the spin degree of freedom is taken into account, a spin coherence develops spon- 
taneously and turns the QH system into a QH ferromagnet, provided the Zeeman effect is 
reasonablly small. Quasiparticles are skyrmions [JTOj , whose existence has been established 
experimentally ||11|| . On the other hand, a pseudospin (interlayer) coherence develops sponta- 
neously in certain bilayer QH systems [12|, [13], |14]] . Some of its characteristic behaviors have 



already been observed experimentally [15, 16|. In a coherent state the phase and number 



differences are observables simultaneously. The phase difference is controlled by applying 
the parallel magnetic field |L5] while the number difference is controlled by applying bias 
voltages to the bilayer system |16| . 



In this paper we analyze skyrmion excitations in monolayer and bilayer QH ferromagnets. 
We use the improved composite-boson (CB) theory M, proposed based on a suggestion due 
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to Girvin 0, Read [§J and Rajaraman et al. 0. The advantages of the improved scheme 
are that the field operator describes solely the physical degrees of freedom representing 
the deviation from the ground state and that the semiclassical property of excitations is 
determined directly by the microscopic wave function. 

In bilayer QH systems we have two types of skyrmions; the spin-skyrmion associated with 
the spin coherence and the pseudospin-skyrmion associated with the pseudospin coherence. 
One skyrmion consists of a pair of charged excitations on the two layers, and hence acquires a 
capacitive charging energy, whichever type of skyrmion it may be. We have predicted |I2|, [TB| 
that the interlayer-coherent state continues to exist even if the electron densities are made 
arbitrarily unbalanced between the two quantum wells, as has been confirmed experimentally 
16| . Furthermore it has been revealed |16] that the activation energy increases as the density 



imbalance increases. We explain this characteristic behavior by evaluating the excitation 
energy of one skyrmion as a function of the density imbalance: The major energy is the 
capacitive charging energy we mentioned. Throughout the paper we use the natural units 
h = c=l. 



2 Bosonization 

We summarize the idea of the improved CB theory applied to the monolayer spin-frozen 
QH state. We denote the electron field by ip(x) and its position by the complex coordinate, 
which we normalize as z = (x + iy)/2lB with £g the magnetic length. Any state |6) at the 
filling factor v — 1/m (m odd) is represented by the wave function, 

6[x] = (0|V>(sbi) • ■■^{x N )\e) = co[z}& LN [x], (1) 

where ©ln[*] is the Laughlin function, 

N 

©lnW =n(^-^) m exp[-^|z r | 2 ], (2) 

r<s r=l 

and u[z] = w(zi, Z2, ■ ■ ■ ,%n) is an analytic function symmetric in all N variables. The map- 
ping from the fermionic wave function &[x] to the bosonic function u[z] defines a bosoniza- 
tion. We call the underlying boson the dressed composite boson and denote its field operator 
by ip(x). The field operator turns out to be the one considered first by Read || and revived 
recently by Rajaraman et al. |J. It follows that 

& v [x] = (OMxt) ■ ■■cp{x N )\e) = co[z). (3) 



The Laughlin state is represented by ©^[aj] = 1. A typical quasiparticle (vortex) state is 
described by ©^[aj] = z ri leading to (f(x)) = z in the semiclassical approximation. This 
is a highly nontrivial constraint, which turns out to determine all semiclassical properties 
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of the vortex excitation. Consequently, the field operator ip(x) describes solely the physical 
degrees of freedom representing the deviation from the ground state, and the semiclassical 
property of excitations is determined directly by the microscopic wave function (p]). 

We now construct the improved CB theory explicitly. We start with the kinetic Hamil- 
tonian for planar electrons in external magnetic field (0, 0, —B), 

where Pj = —idj + eA^ xt is the covariant momentum with A ext = \ejkXkB. Here, £12 = 
—£ 2 i = 1 and £n = £ 2 2 = 0. We denote the electron density by p(x) = ^(x)i[)(x). The 
Coulomb term is 

H c = \ J d 2 xd 2 xV(x - y)g{x)g(y), (5) 

with V{x) = (e 2 /£)|a?| _1 , where g(x) = p{x) — po is the density deviation from its average 
p = (p(x)). It is normalized to vanish, (He) = 0, on the homogeneous ground state. 

At sufficiently low temperature all relevant excitations are those confined to the lowest 
Landau level. The condition is that the kinetic energy is quenched on the state, 

(P X +lPy)^X)\&) =0. (6) 

We call it the lowest-Landau-level (LLL) condition. 

We define the bare CB field by way of an operator phase transformation of the electron 
field ip(x), 

<P(x) = e- l@{x ^(x). (7) 

The phase field Q(x) is chosen to attach m units of Dirac flux quanta 2n/e to each electron 
via the relation, 

EijdidjO(x) = 2-nmp(x). (8) 

When m is odd, <p(x) is proved to be a bosonic operator. The bare composite boson is the 
one familiar in literatures 0, [12], |13"|| . 

We proceed to define the dressed CB field (f(x), 

<p(x) = e~ A{x) (j)(x). (9) 

Here, the hermitian field A(x) is to be determined so that the basic formula (|j) is obtained for 
the wave function. Substituting (||) and (0) into @, the kinetic Hamiltonian is transformed 
into 

Hk = 1 ImS ^V^^-^l^ + ^X^), (10) 
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where <p*(x) = ^{x)e 2A<yX \ with which p{x) = ip^(x)ip(x) = Lp$(x)Lp(x). The covariant 
momentum is 



Vj = -idj + eAf\x) + dj&(x) - idjA(x). (11) 
By a judicious choice of A(x) we are able to bring the LLL condition (§) into a simple form, 

(v x + iv y Mx)\e) = -±Ap(x)\6) = o. (12) 

Indeed, when we choose diA(x) = e^eA^ix) + djQ(x)\, the covariant momentum reads 

Pi = - i9 i ~ ( £ jk + i5 jk )d k A(x), (13) 

from which the LLL condition (|T2|) follows trivially. By using ([8]), the above definition of 
A(x) leads to 

V 2 A(x) = 2nm ( p(x) - — — J . (14) 



2-nm r 

It is easy to see that bare composite bosons feel the effective magnetic field B c s = e~ 1 'V 2 A(x) 
and that B e g vanishes on the homogeneous ground state where (p(x)) = po- It follows that 
the homogeneous ground state is realized only at the filling factor v = 2irp /eB = 1/m, 
where composite bosons undergo bose condensation. Eq . fll4]) is solved as 

A{x)=m I d 2 y\J l -^-^-)g(y), (15) 



2£ B 

in terms of the density deviation g(x) at v = 1/m. It is interpreted that the new CB field 
(HD is obtained by dressing the bare field 4>(x) with a cloud of the effective magnetic field 
generated by A(x), and hence we have called it the dressed field. 

Solving the LLL condition (|P2"D we find that the iV-body wave function ©^[cc] is an 
analytic function as in (R). It is an easy exercise to derive the relation ||, 



■<pHx N )\0) = ©lnM^Ozi) • • -^(x N )\0). 



Because of this relation the function ou[z] in the wave function ([!]) is given precisely by the 
formula (§). 

One might question the hermiticity of the theory ||, since the covariant momentum (|13"D 
has an unusual expression. Analyzing the Lagrangian density we find that the canonical 
conjugate of tp(x) is not i^(x) but i^(x) = iip^ (x)e 2 " A<yX \ It implies that the hermiticity is 
defined together with the measure e 2A ^ x \ Such a measure has arisen since the transformation 
(|9p is not unitary. The covariant momentum ( Jl3| ) is hermitian together with this measure. 
It is instructive to rewrite the kinetic Hamiltonian (pTOT) as 



Hk = | / d\ (^(*)) 'e 2 ^^(x), (16) 



which is manifestly hermitian. 
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3 Monolayer QH Ferromagnets 



We analyze the QH system with the spin degree of freedom. The electron field ip a (x) 
is labeled by the spin index a =f, |. Bare and dressed CB fields <fi a (x) and (p a (x) are 
defined by (0) and (Q), where Q(x) and A(x) are defined by (|8]) and (|1^) with p(x) = 
^ "I" (x)^l(x) = & (x)e 2A ( x }$(x). Here, ^(x) and $(#) are two-component fields made of 
ip a (x) and (p a (x). We decompose the bare CB field into the U(l) field <p(x) and the SU(2) 
field n a (x), <ft a (x) = <j)(x)n a (x), by requiring ^ Q n a ^n a = 1. Thus, 

<p a {x) = e- A{x) <f) a {x) = e- A{x) <j){x)n a {x). (17) 



The field n a (x) is the complex-projective field [fL8|| , whose overall phase has been removed 
and given to the U(l) field 4>(x). 
The kinetic Hamiltonian reads 



Rk = ^m\ d 2 x&{x){V x -iV y ){V x + iV y mx). 



'II 



The LLL condition is 

(V x + iP v )*{x)\6) = $(x)|6> = 0. (19) 

Therefore the CB wave function & v [z] is symmetric and analytic in N coordinates Zi, with 
which the electron wave function is given by &[x] = ©^[^©ln^]- 

We use the index a = x,y,z for the spin component. The spin density is 

S a (x) = ^(x)r a V(x) = ^p{x)s a (x), (20) 

where r a is the Pauli matrix and s a (x) is the normalized spin field, or the nonlinear sigma 
field, 

s a {x) = n\x)T a n{x), (21) 
with n(x) the two-component field made of n^(x) and n^(x). The Zeeman term is 

H z = - l -9*^ B B J d 2 xp(x)s z (x), (22) 

where g* is the gyromagnetic factor and pb the Bohr magneton. Each Landau level contains 
two energy levels for spin- up and spin-down states with the one-particle gap energy g*psB. 
The ground state \go) is unique, whose wave function is 

( \ spin 
J ©lnW, (23) 
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where the two-component spinor is common to all electrons, representing the spin-up polar- 
ization. 

The effective Hamiltonian governing perturbative fluctuations of the sigma field has been 



derived previously [O, W4 



H eS =^ Ps j d 2 x{[d k s(x)] 2 -^g*ti B Bs*(x)y (24) 

The first term represents the spin stiffness with p s = ue 2 / (16\/2ti€£b) ■ Perturbative exci- 
tations are charge neutral. We consider the vanishing limit of the Zeeman term (g* = 0). 
In this case the ground state is given by an arbitrary constant sigma field, s ( x )= const ant. 
All spins are polarized into one arbitrary direction. There exists a degeneracy in the ground 
states. The choice of a ground state implies a spontaneous magnetization, or a QH ferro- 
magnetism. When a continuous symmetry is spontaneously broken, there arises a gapless 
mode known as the Goldstone mode. Quantum coherence develops spontaneously. Actually, 
due to the Zeeman term all spins are polarized into the z axis. As far as the Zeeman effect 
is small enough, the system is still considered as a QH ferromagnet with a finite coherent 
length. 

4 Topological Excitations 

We analyze topological excitations on the QH ferromagnet. The semi classical analysis is 
powerful when the CB wave function &ip[z] is factorizable, 

(| / \ \ spin 
l$)) r (25) 

In this case the one-point function is given by (ip a (x)) = u a (z). Based on the formula (|17D 
it is parametrized as 

e-AWjxWy/pa + e {x)n a {x) = u a (z), (26) 

since |0(a;)| 2 = po + q{ x )- Here and hereafter all fields are classical fields. The Cauchy- 
Riemann equation for (|26| ) yields [ |T3|j , 

q(x) = -vQ (x) + f V 2 lnf 1 + , (27) 

4tt V Po J 

where use was made of (|14|), and Qo(x) is the topological charge density whose explicit form 
is given later. A nontrivial density deviation q(x) is induced around a topological soliton 
according to this equation, which we name the soliton equation. Since the soliton equation 



is a direct consequence of the semiclassical LLL condition (26|), it is interpreted that the 
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density deviation occurs in order to confine topological excitations within the lowest Landau 
level. 

The excitation carries a quantized electric charge, 

q e = — e J d 2 xg{x) = ev J d 2 xQ (x) = uqe, (28) 

as follows from (^7|), where q = f d 2 xQo(x) is the topological charge. We explicitly consider 
solitons with q — 1. One soliton carries a fractional charge, say, q e = |e at v = 1/3, as has 



been experimentally observed [19]. We now show that there are topological solitons, vortices 



and skyrmions, associated with the U(l) and SU(2) parts of the CB field &(x). 



4.1 Vortex Excitations 

The U(l) excitation is generated on the ground state when dkx{ x ) 7^ and n a (x) = constant 
in (|2"E|). The complex-projective field is trivial, n* = 1 and = 0, representing a spin-up 
polarized ground state (^3|). It follows from (p6|) that the topological charge density is 

Q (x) = Qo(x) = —ejkdjdkxix), (29) 

Z7T 

as represents the vorticity. 

The simplest excitation is given by one vortex sitting at x — 0, whose wave function is 
(p5|) with u^(z) = y/poz and ou^(z) = 0. The topological charge is concentrated at the vortex 
center, Q v (x) = 5(x). An approximate solution of the soliton equation (^7]) reads 

^-A i+ t-iy^" B - « 3o) 

A vortex is a hole made in the condensate of composite bosons. 



4.2 Skyrmion Excitations 

The SU(2) excitation is generated on the ground state when x{ x ) = constant and dkn a {x) ^ 
in (|26|). The complex-projective field is solved as 

n a (x) = l £M (31) 

with arbitrary analytic functions u a (z). It is known to describe skyrmions [18|]. We call them 
spin-skyrmions since they are associated with the spin coherence. The simplest excitation is 
given by one skyrmion with scale k sitting at x — 0, whose wave function is 

(\ spin 
l) 2 ) ©ln[4 (32) 
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The skyrmion is reduced to the vortex in the limit k — > 0. The scale k is to be fixed 
dynamically. The normalized spin fl2~T|) is calculated from this wave function as 

- 2 - ( w 

r 2 + (£ b k) 2 ' 



s x sky = Vl-(s z ) 2 cose, 3 v Ay = -y/l-(a*)*Bin9, s s \ y = l , ), B *{ 2 . (33) 



The spin flips at the skyrmion center, s = (0,0, —1) at r = 0, while the spin-up-polarized 
ground state is approached away from it, s = (0, 0, 1) for r ^> k£b- 

The soliton equation (|2?|) follows together with the topological charge density, 



Q (x) = Q?(x) = -±.J2^kd j (n a *d k n a ). (34) 



The topological charge is shown |L8[ to be identical with the Pontryagin number whose 
current density is 

Q p Jx) = ^- £ahc e^ x s a d»s b d x s c . (35) 



Approximate solutions of the soliton equation (pTj ) are constructed in the two limits, the 
large skyrmion limit (k ^> 1) and the small skyrmion limit (k <C 1). First, in the large limit 
we can solve (|27|) iteratively, where the first order term is 

Q*y(x) ^ -uQ p (x) = - U - [r ^ B)r (36) 

This agrees with the formula due to Sondhi et al. [ 10| . However, in the small limit the 
topological charge Qq(x) is localized within the core. Indeed, we have Qq(x) — > 5(x) as 
k, — > in (|36|), with which the solution of the soliton equation ( f27j ) is approximated by 
the vortex configuration g vor (x) in (0). This is what we have expected since the skyrmion 
wave function is reduced to the vortex wave function in the limit k — > 0, where there is no 
distinction between the U(l) and SU(2) excitations. Since the vortex may be considered as 
a small skyrmion limit, we do not make a clear distinction between them in what follows. 

5 Bilayer QH Systems 

We proceed to analyze bilayer QH systems. In this section we freeze the spin degree of 
freedom. We are interested in the so-called Halperin (m, m, m) phase ||20|| , or the interlayer- 



coherent phase, in which an interlayer quantum coherence develops spontaneously |l2[ \L3 
We denote the electron field at the layer a{= 1, 2) by if) a . We call the layer a = 1 the front 
layer and the layer a = 2 the back layer. It is convenient to introduce the pseudospin SU(2) 
structure by considering a two-component electron field ty(x) made of ipi and ip2- We use the 
pseudospin index a = 1, 2, 3 and define the pseudospin density and the normalized pseudospin 
field by ( |2"0"D and ([2"1~D, respectively. In the interlayer-coherent phase, the dressed CB field 
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$(x) is defined precisely as in the monolayer QH ferromagnet. The kinetic Hamiltonian of 
composite bosons is formally identical to that of the monolayer QH ferromagnet, and given 
by fll8|). Consequently, the LLL condition is identical to (|19|), and the Hilbert space made of 



the states in the lowest Landau level is identical to that of the monolayer QH ferromagnet. 
We make a full use of this mathematical identity (13| to investigate the interlayer-coherent 
phase. 



The Coulomb interaction is summarized into the two terms [ 13 ], 

H c = \f d 2 xd 2 xV + (x - y) e (x)e(y), (37) 
H c =2 J d 2 xd 2 xV'{x - y)AS 3 (x)AS 3 (y), (38) 

where V ± (x) = (e 2 /2s)(\x\~ 1 ± (|a;| 2 + d 2 ) -1 / 2 ) with the interlayer separation d: g(x) is 
the deviation of the total electron density. The Coulomb energy H@ is the driving force to 
realize the QH system. The term describes the capacitive charging energy between the 
two layers. The deviation of the pseudospin density from the ground-state value is 

AS a (x) = ^[p(x)s a (x)-p s^, (39) 

where s a (x) is the normalized pseudospin and Sq is its average value in the ground state: 
See (f43|). On the other hand, the tunneling term is given by 

H T = -A SAS [ d 2 xAS\x). (40) 



The tunneling term produces an energy gap Asas between the symmetric and antisymmetric 
states. The pseudospin (interlayer) coherence develops provided the capacitance term 
and the tunneling term Ht are reasonably small. 

There is one additional degree of freedom in the bilayer system. By applying bias voltages 
to the two layers, we can freely control the electron density p$ in each layer |16| . Namely, the 



direction of the pseudospin polarization is controllable [fL3|| . Accordingly, the ground state 
wave function is given by 

^[^^(^l)" 6 ^], (41) 

where the pseudospinor is common to all electrons and the parameter do is a real constant 
(|c"o| < 1); the index ppin denotes the pseudospin. Note that the ground state (|4l]) has been 
chosen so as to minimize the tunneling energy. The density ratio between the two layers is 

4 = ^ («) 
Po 1 - ff o 
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The normalized pseudospin is calculated from the wave function ( IIP as 



sl(x) = yjl- crl s 2 (x) = 0, s 3 (x) = a . (43) 

The balanced configuration is realized when o"o = 0, where the symmetric state is the ground 
state and described by (|4"lD with <Jo = 0. 

We analyze charged excitations on the ground state (|4lD . The easiest way is to use 
a mapping between the bilayer pseudospin state and the monolayer spin state jl3| . The 
mapping is established by the matrix, 

T - — ( V1 + a ° V1 ~ a ° \ (44) 

It is easy to see that T a & ppm = ©g Pin , where & ppm is the ground state ([y]) in the unbalanced 
bilayer system and <S spm is the ground state (^) in the monolayer QH ferromagnet. All 
excitations are mapped by this transformation between the two systems. Consequently, the 
skyrmion excitation on the ground state (|4~ID is obtained from that on (|23"D as = T^& S J^ 
with (0), or 

This is the pseudospin-skyrmion associated with the pseudospin coherence. The normalized 
pseudospin ( pi]) is easily calculable from this wave function with the aid of (|3T|). It turns 
out that one skyrmion excitation consists of two parts, one on the front layer and the other 
on the back layer. 



6 Bilayer QH Ferromagnets 

We finally analyze the bilayer QH system with spins included. The lowest Landau level 
contains four energy levels corresponding to the two layers and the two spin states. In 
the balanced configuration the ground state is the symmetric and spin-up state, which is 
separated by the Zeeman gap g*\iB^ and/or the tunneling gap Asas from the one-particle 
excited states. In the unbalanced configuration at v — 1/m the QH state is given by 




as a tensor product of the spinor and the pseudospinor. When both of the interactions are 
much smaller than the Coulomb energy, two kinds of quantum coherences develop sponta- 
neously upon this ground state. 

We consider the filling factor v = 1 and v = 2. There are two cases. (A) When the 
Zeeman gap is smaller than the tunneling gap, spin-skyrmions are excited at v = 1 while 
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pseudospin-skyrmions are excited at v = 2. (B) When the Zeeman gap is larger than the 
tunneling gap, pseudospin-skyrmions are excited at v — 1 while spin-skyrmions are excited 
at v = 2. 

We are interested in the activation energy A act (<7o) as a function of the density imbalance 
parameter a . Though A act (<7 ) is sensitive to impurities in samples, the difference A act (cr ) — 
A act (0) is not so. The excitation energy of the spin-skyrmion depends on <jq only through the 
capacitance term (|3"8"D, but that of the pseudospin-skyrmion depends also on the tunneling 
term fl4"0"|). The capacitive charging energy arises since one skyrmion consists of two charge 
excitations, one in the front layer (q\) and the other on the back layer {(£,). It carries the 
total electric charge ( |2%D with q e = q\ + q\ = ve. 

We analyze spin-skyrmions. Though the normalized pseudospin ( |43|) is not affected by 
the spin-skyrmion excitation, it induces a deviation of the pseudospin fl39f ), 



<AS 3 (*)> = \{q^{x))s% = ^g sky (x). (47) 

The charge difference is 

q\ - q 2 e = _2e / d 2 x(AS 3 {x)} = ua e, (48) 



independently of the details of excitations. Hence, the charge on the front layer is q\ = 
\ve(X + o'o)) while that on the back layer is q 2 e = \ve{\ — <j ). The ratio q\jq\ is the same 
as the ratio (^) of the electron densities on the two layers, as is expected. The capacitive 
charging energy (|38|) is estimated as 



eV 2 



(He) = «c^, (49) 

with a numerical constant ac at a fixed value of djls (the ratio of the interlayer distance 
d and the magnetic length 1%). Its value depends on the details of the excitation. The 
pseudospin-skyrmion is similarly analyzed, and the same formula as (^) is obtained for its 
charging energy. 



We have compared the above results with the experimental data taken from Ref.|L6 
where g* \ibB / Asas — 1/4 at B = 5 Tesla. See Fig.^. At v = 1 we have fitted the 
activation-energy data by assuming a spin-skyrmion excitation with size k = 2, where the 
relevant parameters are d = 231 A and £b = 120 A. At v = 2 we have fitted the data by 
assuming a pseudospin-skyrmion excitation with size k = 1.65 and charge 2e, where the 
relevant parameters are d = 23lA, £g = 228A and Asas = 6K. The agreements are quite 
good, as would imply that the present picture of the bilayer QH ferromagnet is basically 
correct. Detailed analysis will be published in a separate paper |2T] . 
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Figure 1: The excitation energy (in K) of one skyrmion is calculated as a function of the 



imbalance parameter ctq, and compared with the experimental data taken from Ref.|T^]. We 
have adjusted the value at the balanced point (ctq = 0) to the experimental one. 



7 Discussions 



We have studied both the monolayer and bilayer QH ferromagnets based on the improved 
CB theory We have investigated soliton excitations confined to the lowest Landau level. The 
semiclassical approximation is powerful when the CB wave function & v is factorizable. This 
allows us to analyze quasiholes successfully. It is quite difficult to make a similar analysis 
of quasielectrons, for which (5^ is not factorizable. In comparing our theoretical results 
with experimental data, we have made a physical assumption that the excitation energy of 
one quasielectron is approximately the same as that of one quasihole. In order to discuss 
quasielectrons, it would be necessary to make the LLL projection after exciting them. We 
hope to discuss on this point in a future paper. 
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Ministry of Education, Science, Sports and Culture. 
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